Towards a black body for near-field thermal radiation 
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We study the near-field heat exchange between hyperbohc materials and demonstrate that these 
media are able to support broadband frustrated modes which transport heat by photon tunnelling 
with a high efficiency close to the theoretical limit. We predict that hyperbolic materials can be 
designed to be perfect thermal emitters at nanoscale and derive the near-field analog of the black- 
body limit. 
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A black body is usually defined by its property of hav- 
ing a maximum absorptivity and therefore also a maxi- 
mum emissivity by virtue of Kirchhoff's law The en- 
ergy transmission between two black bodies having differ- 
ent temperatures obey the well-known Stefan-Boltzmann 
law. This law sets an upper limit for the power which 
can be transmitted by real materials, but it is itself a 
limit for the far-field only, since it takes only propagating 
modes into account. In terms of the energy transmission 
between two bodies the black body case corresponds to 
maximum transmission for all allowed frequencies uj and 
all wave vectors smaller than uj /c, where c is the vacuum 
light velocity. This means that all the propagating modes 
are perfectly transmitted across the separation gap. 

In the near- field regime, i.e., for distances smaller 
than the thermal wavelength Ath = hc/k^T (27r/i is 
Planck's constant, /cb is Boltzmann's constant, and T 
is the temperature) the radiative heat flux is not due to 
the propagating modes, but it is dominated by evanes- 
cent waves I2tJ1 and especially surface polaritons as 
confirmed by recent experiments [5-ll|. The common 
paradigm is that the largest heat flux can be achieved 
when the materials support surface polaritons which will 
give a resonant energy transfer restricted to a small fre- 
quency band around the surface mode resonance fre- 
quency 0, 3, 13 1- Many researchers have tried to 



find materials enhancing the nanoscale heat flux due to 
the contribution of the coupled surface modes by using 
layered materials lj,ll5|, doped silicon Imll7[, metama- 
terials Il8l-l20l . phase-change materials |2l[ and recently 
graphene [22|. 

In the present work the aim is twofold: (i) We show, 
that materials supporting a broad band of evanescent 
frustrated modes can outperform the heat flux due to 
surface modes. This provides new possibilies for design- 
ing materials giving large nanoscale heat fluxes which 
could be used for thermal management at the nanoscale 
for instance, (ii) We derive a general limit for the heat 
flux carried by the frustrated modes and show that it 



is, in fact, the near-field analog of the usual black body 
limit. For the evanescent modes a near field analog of 
a black body can be defined in the sense that the en- 
ergy transmission coefficient must be equal to one for all 
frequencies and all wave vectors larger than w/c. 

With today's nanofabrication techniques it is possi- 
ble to manufacture artificial materials such as photonic 
band gap materials and metamaterials which exhibit very 
unusual material properties like negative refraction [23|. 
Due to such properties the y a re considered as good can- 
didates for perfect lensing j24L [2H | , for repulsive Casimir 
forces |26l - l29| and enhanced or tunable radiative heat fiux 
at the nanoscale to mention a few. 

There exists a class of uniaxial metamaterials for which 
the permittivity and permeability tensor elements are not 
all of the same sign [33|. In particular, for such materi- 
als the dispersion relation for the solutions of Helmholtz's 
equation inside the material is not an ellipsoid as for nor- 
mal uniaxial materials [s^l but a hyperboloid [s^. For 
this reason such materials are also called hyperbolic ma- 
terials. These materials have already been considered for 
super-resolution imaging [s^ and enhanced thermal con- 
ductivity inside the material itself [s^l for instance. Here 
we focus on the heat fiux between two bodies consisting 
of hyperbolic materials showing that these materials can 
be used to realize a black body at the nanoscale. 

The nanoscale heat flux <I> between two semi-infinite 
materials separated by a vacuum gap of thickness d 
having the temperatures T 4- AT and T can be writ- 
ten in terms of the quantum of thermal conductance 
-K^k^T/Zh [s^ and the sum over all transversal modes 
as bol 
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where k = {kx,kyY is the lateral wave vector. The mean 
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transmission coefficient introduced here is defined as 

' /duu2e«(e"-l)-2 ^ ' 

where u = hw/k^T is a rescalcd frequency or energy 
and Tj(Lu,K) is the energy transmission coefficient for s- 
and p-polarized modes in (w, n) space. For isotropic and 
anisotropic semi-infinite materials the cxphcit expression 
of this energy transmission coefficient is known and can 
be found in Refs. d H El j . 

From Eq. ([T]) it is obvious that the heat flux $ can be 
increased on the one hand by increasing the mean trans- 
mission coefficient for a transversal mode with a given k 
or by increasing the number of modes contributing sig- 
nificantly to the heat fiux on the other. As discussed 
in Refs. [3, ll^i the number of propagating modes con- 
tributing to the heat flux is limited due to the fact that 
the dispersion relation of these modes -I- k^Q = ixP' jc? 
restricts for a fixed frequency the modes to a circle 
K < w/c. For dielectric materials the frustrated inter- 
nal refiection modes [42j which can also contribute to the 
heat flux for distances d < Ath the dispersion relation 
reads + kli = tuP' jt? so that in this case the con- 
tributing modes are for a fixed frequency restricted to a 
circle k < -y/ew/c. 

in presence of surface resonant modes such as surface 
phonon polaritons 

It is known that in presence of surface resonant modes 
such as surface phonon polaritons 4^ 44 1 in the infrared 
the heat flux at nanoscale can be increased by several 
orders of magnitude (see Fig.[T]). In terms of the 

mean transmission coefficient Tj it can be shown (soj 
that the transmission of these modes is very small (less 
than one percent, see Fig. 3 of Ref. [l^), but the con- 
tributing modes are for very small distances restricted to 
K < log[2/Im(e)]/2(i. First, this means that the number 
of contributing modes goes like 1/(P for very small dis- 
tances so at very small distances the surface phonon po- 
lariton contribution to the heat flux becomes larger than 
the contribution of the frustrated total internal reflec- 
tion modes. Therefore the surface phonon polaritons pro- 
vide the dominant heat flux channel at the nanoscale de- 
spite the fact that their transmission coefficient is rather 
small (39I I . Second, the upper limit of the wave vectors k 
for th e contr ibuting modes also depends on the intrinsic 
losses 



45[ of the material. Note, that there is a funda- 
mental limit restricting k to k < tt/o where a is for polar 
materials given by the lattice constant so that there is 
also a definitive limit for the nanoscale heat flux [46. 47|. 

Now we are in a position to discuss the heat flux for 
hyperbolic materials. Such materials are a special class 
of anisotropic uniaxial materials. The dispersion rela- 
tion for extraordinary waves in uniaxial materials with 
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Figure 1: Heat flux <I> from Eq. ([T]) between two semi- infinite 
SiC bodies over distance normalized to the black body value 
$BB for temperature T = 300 K and AT = 1 K. The plot 
also shows the contribution of the propagating and frustrated 
modes and the surface phonon polaritons (SPhP). 



an optical axis normal to the surface reads 34 1 
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where e± and e|| are the permittivity perpendicular and 
parallel to the surface, respectively, having an opposite 



sign [33|. Hence, the dispersion relation describes not 



an ellipse in the (k, fcj,) plane but a hyperbolic function 
as illustrated in Fig. [Hb). This means that in principle 
without losses there is no upper bound for the k which 
can fulfill the dispersion relation ([3]). In the following 
we will call such modes hyperbolic modes (HM). If we 
consider the nanoscale heat flux for such materials then 
the number of contributing modes is only restricted by 
the intrinsic cutoff in the energy transmission coefficcnt 
which is 7}(w, k) oc exp{—2Kd) for k ^ co/c (see Eq. (O). 
It follows that the heat flux due to the HM scales for 
small distances like l/d^ as the contribution of the surface 
phonon polaritons. When e|| > and e_L < (ey < and 
e± > 0) the HM are propagating inside the hyperbolic 
material for k > w/c (k > e±uj/c). But both kinds of HM 
are evanescent inside the vacuum region, i.e., for small 
distances d <C Ath they are a special kind of frustrated 
internal reflection modes. As for usual frustrated modes 
we can expect a large mean transmission coefficient for 
the HM up to K ~ l/2d if they exist in a broad frequency 
band. Since the mean transmission coefficient for the 
surface phonon polaritons was found to be very small, 
because these modes are restricted to a small frequency 
band around the surface mode frequency, we can expect 
to have a heat flux due to the HM larger than that of the 
coupled surface phonon polaritons. 

In order to see whether our expectation is true, we 
consider a simple example of a uniaxial material which is 
given by periodically positioned nanowires |48| of SiC as 
sketched in Fig. HJa) . In the long wave-length limit [i^ 
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Figure 2: (a) Sketch of a hyperbolic material consisting of 
a periodical array of nanowires. (b) Hyperbolic dispersion 
relation in such materials for a fixed frequency. The solid 
line corresponds to ey < 0, ex > and the dashed line to 
£11 > 0, ex < 0. 



we can describe the effective material properties by the 
Maxwell- Garnett expressions [i^ 
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Here / is the volume filling fraction of the SiC wires and 
e{uj) is the permittivity of SiC As pointed out in 

Rcf. the effective medium theory is limited to dis- 
tances d > o/tt, where a is the lattice constant of the 
nanowires. For smaller distances a nonlocal model for 
the permittivity would be necessary. 

Since the optical axis of the considered material is per- 
pendicular to the surface the s- and p-polarized modes 
decouple and we have only extraordinary waves for p po- 
larization (34| . Then the energy transmission coefficient 
Tj for two identical halfspaces which enters into Eq. ([2]) 
reads 
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for j = {s, p} where Dj : 
like denominator, k^o — \JuP- jc?- 
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introducing kg = t\\LJ^ j (? — and fcp = fc^i = 
e||a;^/c^ — ^k^. From the expression for the reflec- 



tion coefficients it becomes obvious that only the extra- 
ordinary waves represented by rp are sensitive to e_L and 
e||. Hence only the p-polarized part can support HM 
because here we are just considering non magnetic ma- 
terials. For magnetic materials also the s-polarized part 
can support HM. 

In Fig. [3ja) we plot Re(e||) and Rc(e^) in a fre- 
quency range around the transversal and longitudinal 
phonon frequency of SiC (wt = 1.495 • 10^'*rad/s and 
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Figure 3: (a) Real part of ey and ex defined in Eqs. ((4| and 
((S]) with / — 0.1. The two gray areas mark the frequency 
region Ai and A2 where both permittivities have different 
signs, i.e., where the modes fuIlfiU a hyperbolic dispersion 
relation, (b) Transmission coefHcient 7^(a;,K) choosing d = 
100 nm for the same frequency region and the same volume 
filling fraction showing high transmission for regions Ai and 
A2. (c) Heat transfer coefficient for two hyperbolic materials 
with / = 0.1, 0.2, 0.3 normalized to the value for two flat SiC 
half spaces (T = 300 K). Inset: Ai and A2 over /. 



CJL = 1.827 • 10"rad/s) using / = 0.1. The regions Ai 
and A2 where the permittivities have different signs are 
highlighted. Within these frequency bands we have HM. 
In Fig. ^h) it can be seen that the transmission coef- 
ficient Tp{u},K) is close to one for these modes up to 
nd sa 1/2. Furthermore, we observe in Fig. [SJb) that 
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for the chosen parameters we have no coupled surface 
phonon polaritons in our hyperbolic material so that the 
heat flux is only due to the usual frustrated modes and 
the HM. It can be shown that for filling factors / < 1/3 
there are no surface modes in our structure regardless of 
the material used for the nanowires. 

By increasing the filling factor in our structure the two 
regions Ai and A2 supporting HM get broader (see inset 
of Fig. ^c)) resulting in larger heat fluxes in the small 
distance regime around d = 100 nm. It can be seen in 
Fig. [3Jc) that in this distance regime the heat flux due 
to HM can be larger than that by surface phonon polari- 
tons for two SiC media. This makes clear that the goal 
for designing a hyperbolic material with elevated heat 
fluxes is to make the bands for HM as broad as possible. 
This can be done by changing the constituents of the hy- 
perbolic material and by considering different geometries 
like for example different nanowire structures or nanolay- 
ered materials. But there is a general limit for the heat 
fluxes which can be achieved by such materials. 

Now we derive the general theoretical limit for the heat 
flux due to HM. That means we will start with the gen- 
eral expression for the heat flux in Eq. ([l} using the re- 
flection coefficients in Eqs. ([7]) without assuming a spe- 
cial realization of the hyperbolic material as the example 
of the nanowired material. The only assumptions made 
arc that we have a uni-axial local medium with its opti- 
cal axis parallel to the surface normal. Before we start, 
we remind the reader that the usual black body limit 
can be obtained from Eq. ((T|) by considering the limits 
Re(ej_), Re(e||) — >■ 1 and Im(e_L), Im(e||) — ?> for aU fre- 
quencies. Then the Tp — and therefore 7^ — > 1 for 
K < Lo/c and 7^—5-0 for k, > u/c. 

With a similar reasoning we will now derive the the- 
oretical limit for near-field radiation by the HM. For 
convenience, we consider only the case Rc(e||) > 
and Re(e^) < 0. Let us first again assume that 
Im(e_L), Im(e||) — > for all frequencies then we can re- 
cast the refiection coefficient for p-polarized waves as 
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where we have introduced the symbol a = ■\/e|||e_L|7/7± 
with 7^ = - w^/c^ and 7]^ = \e±\uj'^/c'^ + k^. Since 
we consider only the frustrated total internal reflection 
modes {k, > uj/c) a is a real number. Hence, as can 
be expected for total internal reflection modes we have 
|rpp ~ 1. That means one can express the reflection 
coefficient as rp = exjp[iip) with tan(iy9) = Im(rp)/Re(rp). 
Plugging this reflection coefficient into the transmission 
coefficient for the evanescent modes in Eq. ([5]) wc find 
that 7^ = 1 for |fczo|<i ^ 1. On the other hand for 
l^zoM 3> 1 we find that Tp has a maximum for rp = ±i. It 
follows from Eq. ^ that the condition for having perfect 
transmission is = 1 which translates in the quasi-static 
regime (k 3> i^/c) into e|||ei| = 1. 



Finally, we can derive the maximum heat flux due to 
the frustrated HM by plugging rp = ±i into the expres- 
sion for the transmission coefficient in Eq. © and evalu- 
ating Eq. (H)) for the heat flux. Then, we find a maximum 
heat transfer coefficient ($ = hAT) of 
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This expression constitutes the black body limit for the 
near-field radiation by HM, since it is due to perfect 
transmission (7^ = lfora;/c<K<g; l/{2d)) by pho- 
ton tunneling. Note, that this limit is by a factor ln(2)/2 
smaller than the upper limit found in Ref. [soj . which 
was derived by making an assumption on the cutoff in 
the K integration and is therefore ambiguous. 

Let us now compare the upper limit in Eq. ^ with 
known results. For T — 300 K the heat transfer co- 
efficient for a black body is about 6.1Wm~^K~^. Ac- 
cording to Eq. its near-field counterpart (/imax.p ~ 
3.3 • lO^Wm-^K-i) is about 5.4 • 10^ times larger at 
d = 10 nm. For the same distance the heat transfer coef- 
ficient between two gold plates is about 1157Wm^^K^^ 
and between two SiC plates one gets 9200Wm^^K^^. 



Hence h„ 



is by a factor of about 285 or 36 larger. 



rcsp. Finally, we compare the limit for the heat transfer 
coefficient due to photon tunneling of HM from Eq. © 
with a numerically found upper limit for Drudc materi- 
als for the heat transfer coefficient due to the coupling 
of surface plasmon polaritons (see parameter set A in 



Ref. [511) which is about 3 • 10^ times the black body 
value at d = 10 nm. That means the limit in Eq. © 
is by a factor 5/3 larger than the value one can get by 
tuning surface mode resonances. Hence, the design of hy- 
perbolic materials constitutes a new way for an efficient 
enhancement of the nanoscale heat flux. 

In conclusion, we have studied the near-fleld heat 
transfer between hyperbolic materials and we have shown 
with a simple example that the broadband spectrum of 
frustrated modes supported by these media allows us to 
get heat fluxes which are larger than the heat flux due 
to narrow-band coupled surface polariton modes. We 
have derived a theoretical limit of the heat flux due to 
hyperbolic modes and demonstrated that it constitutes 
the near-field analog of the black-body limit. This opens 
up new possibilities to achieve large heat fluxes at the 
nanoscale by the design of hyperbolic materials. 
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